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Formulating a perfect fluid fiUed sphericaUy symmetric metric utilizing the 3+1 formal- 
ism for general relativity, we show that the metric coefficients are completely determined 
^ ■ by the mass-energy distribution, and its time rate of change on an initial spacelike 

' hypersurface. Rather than specifying Schwarzschild coordinates for the exterior of the 

' collapsing region, we let the interior dictate the form of the solution in the exterior, and 

f**^ , thus both regions are found to be written in one coordinate patch. This not only alle- 

viates the need for complicated matching schemes at the interface, but also finds a new 
coordinate system for the Schwarzschild spacetime expressed in generalized Painleve- 
' GuIIstrand coordinates. 
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W)' The traditional approach to the analysis of gravitational collapse follows that de- 

vised by Oppenheimer and Snyder,^ whereby the Einstein field equations are solved 
for the interior, matter-filled region without consideration of the exterior. Whence a 
■ solution to the interior is found, the Israel-Darmois matching conditions are utilized 

to "glue" the interior spacetime to an appropriate exterior spacetime, commonly 
Schwarzschild. Throughout the process the two regions are considered as separate 
entities, mainly as they are described by different coordinate systems. 

We introduce a different approach, whereby the spacetime is established as an 
initial/boundary value problem, with the interface between the two regions of the 
spacetime being a free-boundary. This enables us to describe both regions of the 
spacetime under a single coordinate patch by simply letting the energy-momentum 
variables go to zero at some finite coordinate radius. In this talk we use our formal- 
ism to describe the gravitational collapse of a spherically symmetric perfect fluid. 

As we are setting up an initial value problem, an ideal starting point is the ADM 
system of equations. This will enable us to establish an initial spacelike hypersurface, 
with perfect fluid for r <ra and vacuum for r > rg, where rg is some radius on the 
initial slice. The system can then be evolved forward in time to describe the entire 
collapse process. Furthermore, in order to describe both regions of the spacetime 
utilizing a single coordinate system, one requires that the observer have a flnite 
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radial velocity such that this observer will pass through the interface between the 
two regions. In ADM language, this implies a non-zero, radial component of the 
shift vector, f3(t,r). We therefore begin with an arbitrary, spherically symmetric 
line element expressed as 

dS^ = -a'^dt^ + (1 + 5)"^ {Pdt + drf + r^dn^, (1) 

where Q;(t, r) is the lapse function, £{t,r) > —1 is an arbitrary function which 
reduces to the energy function of the Lemaitre-Tolman (LT) metric, ^'^ and dil^ is 
the metric of the two sphere. 

By putting this line element through the ADM equations (for details see^'^), 
one can derive the reduced field equations which are a coupled system of first order 
differential equations. We define a "mass" function'* 

M := 47r [ p{t, s)s^ds, (2) 

where p(t, r) is the mass-energy density. The lapse function is related to the density 
and pressure through Euler's equation 

drP=-ip + P)dr{\na). (3) 

The solution of this equation requires the specification of an equation of state (EoS) 
which relates the density to the pressure, P. Thus, given an EoS, the system reduces 
to the line element along with two equations. 



dS^ = (1 + ^^)' 



n / 2M\ o 2M n o n 

-a^ I 1 ]dt^ + 2a\l \- £ dtdr + dr^ + r^dfl^ 

r 



(4) 



Cn£= '^^^^^\ —+£ drP and £^M = 47rPr\ —+£. (5) 
p + P \ r \ r 

Here, denotes the Lie derivative with respect to the unit normal vector which, 
when acting on a scalar, takes the form £„'0 = {dtip — /3drtp). Once an EoS 
is specified, equation (O is solved and hence the lapse function can be written in 
terms of the density, and thus the mass. Therefore, equations (O are two equations 
for two unknown functions M and £. 

As per our aim, these equations describe both the interior perfect fluid region 
and the exterior vacuum region in the one coordinate patch. To see this we simply 
let the pressure and density vanish external to some finite radius. Equation |(2J) 
then implies that the mass function is constant, which further implies the right 
hand equation in JSJ is trivially satisfied. Equation implies the lapse function is 
simply a function of the temporal coordinate, and utilizing coordinate freedom the 
lapse can be set to unity without loss of generality. The resulting system is what 
we call the generalized Painlevc-GuUstrand (GPG) line element as the special case 
{£ = 0) is the Painlcve-GuUstrand line element. 



''We note in the dust limit, M becomes the famiUar mass of the LT solution, and in the vacuum 
limit becomes the Schwarzschild mass. 
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The GPG class of solutions comprise a family of coordinate systems for the 
Schwarzschild spacetime. The coordinate transformation between this class of so- 
lutions and Schwarzschild coordinates, (i,r,0,ip^, is given by the solution to the 
coupled differential equations 



{dity = l + £ and 1 ]drt^\ +5, (6) 



where t = t{t,r). One can show a solution exists to these equations, and thus the 
coordinate transformation is always valid by utilizing the integrability conditions, 
which are satisfied providing the left hand equation in lO is satisfied (for details 
see''). 

Reverting back to the full system of equations in the interior with matter (|2I5|I . 
we can transfer these into diagonal coordinates, {T,R,d,(j)), such that the general- 
ization from the LT metric for dust becomes obvious. By letting r — r{T, R) such 
that 

{drrf = a^{^+£\, (7) 



the line element and equations ^ become 



n2 



rf52 = -a^dT^ + ^^^dR^ + rdn\ (8) 
1 ^ c 



1 I p /2M 2M 

(dnr) (dr£) = —^aJ + £ dnP and OrM = 4TrPr^a\ +£, (9) 

p + P \ r \ r 

and equations |2Jl and Q are suitably dealt with. The reduction to the LT dust 
models is clear, again using the coordinate freedom that the lapse function becomes 
a function of time, and can thus be set to unity without loss of generality. 
A number of extensions of this work are currently under investigation: 



Searching for exact solutions of equations ©-lEl- 

The analysis and determination of shell-crossing singularities which exhibit 

themselves as fluid shock waves in this coordinate system. 

The relaxation of the perfect fluid condition to allow for more realistic 

matter sources, enabling the study of diffusion processes and anisotropic 

stresses. 

A relaxation of the symmetries of the geometry to allow for quasi-spherical 
symmetry or axial-symmetry. 
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